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Abstract 

A shape of lateral distribution for charged particles in events with energy above 10 19 eV is considered. 
Two methods were used for individual LDF parametrization. In the first approach, the index of power was 
determined for generalized Greisen-Linsley approximation. In second, mean square radius of the shower 
was determined for approximation proposed by Lagutin et al. Comparison of resulted parameters is pre- 
sented for individual events arrived from different celestial regions — Galactic planes and the region with 
increased flux of particles with Eq > 10 19 eV (according to Yakutsk array): 1.7h — 3.7h right ascension ; 
45° - 60° declination. 



L INTRODUCTION 

The knowledge of the lateral distribution 
function (LDF) of charged particles from ex- 
tensive air shower (EAS) is vital for experi- 
ments in the field of ultra-high energy cosmic 
ray (UHECR) studying. It is LDF that defines 
main shower parameters such as p 60 o (charged 
particle density at the distance 600 m from the 
core) and thus — primary energy. 

In this paper we consider parameters of indi- 
vidual LDFs resulted from revision of high en- 
ergy events registered at the Yakutsk EAS array. 
The aim of this work is to trace possible corre- 
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lation between parameters of individual showers 
and their arrival directions on the sky, especially 
for Galactic planes and for the region with sig- 
nificantly increased UHECR flux, detected by 
Yakutsk group Q]|. 

H. ESTIMATION OF LATERAL DISTRI- 
BUTION PARAMETERS FOR INDIVIDUAL 
SHOWERS 

For the analysis we selected showers with 
E > 10 19 eV, with zenith angles 9 < 60° and 
with core lying well within the boundaries of the 
array, to make sure that shower core is found 
correctly. 

At the Yakutsk EAS array, approximation 
proposed by Greisen [0 is used for primary 



data processing: 

p(r) = M 



r 

Ro 



-1 



1 + 



Re 



(b)+l 



, (1) 



where R is Moiere radius and slope parameter 
(b) = -1.38 - 2.16 • cos9 - 0.15 • lgp 600 . 

In the work by Glushkov et al 0, an updated 
approximation was proposed, that demonstrated 
better description of experimental points at large 
distances from the core (r > 1000 m): 
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Figure 1: Ab vs Galactic latitude. Marked stripe 
±10° is Galaxy plane. 



x 



2000 



-3.5 



(2) 

where (b) = 2.6 • (1 - cosfl) - 3.242. 

In equations ([T]) and Q the slope parameter 
(b) is derived from average LDF. While it de- 
scribes most of showers quite well, it certainly 
fails doing so in dozen number of events. Dur- 
ing revision we performed x 2 ~fitti n g of func- 
tions ([T} and ([2]) normalized to p 60 o on experi- 
mental data for each selected shower with free 
parameters p 60 o and b. 

The value Ab = \{b) — b\ could give a hint 
of possible astrophysical aspect of the slope 
parameter in functions ([T]) and Q. As seen 
on FigjTJ comparison to Galactic coordinates 
showed no correlation between Ab and Galaxy 
plane. 

To trace possible dependency on astrophys- 
ical coordinates, we selected a "stripe" of 15° 
width along declination and divided it into 
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Figure 2: Ab vs astrophysical coordinates. 

"chunks" of 2h each along right ascension. Such 
a stripe was selected to exclude zenith-angular 
dependency. Averaged Ab values in each chunk 
are presented in table [I] in comparison to aver- 
aged value in the rest chunks of the stripe. 

A. Scaling approach 

A one-parametric scaling representation of 
charged particle lateral distribution was pro- 
posed by Lagutin et al : 
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Table I: Averaged Ab values in the region of interest (see Fig. [2]) 
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here i? ms is mean square radius of electrons. 
This function was obtained with respect to 
nuclear cascade process in the shower JH. 
Since the main classification parameter for the 
Yakutsk array is p 60 o, we used ([3]) normalized 
to p 6 oo- We calculated i? m . s . for each shower in 
our selection using ^-minimization. On Fig. [3] 



there are shown R m s values obtained for indi- 
vidual events compared to zenith angle. It is 
clear, that these values significantly exceed pre- 
dicted in the work [4], though one can note dis- 
tinct zenith- angular dependence. 

We constructed average LDFs for three 
zenith- angular intervals: — 30°, 30 — 45° and 
45 — 60°. Results can be found in table [TT] and 
Figj4j It is seen from the table, that resulted 
R m s . values contradict to theoretical predictions 
from the work JH. 

Obtained i? ms values did not allow us to 
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Figure 3: R m . s . of electrons in individual showers 
compared to sec 9. 

Table II: Parameters for average LDF obtained for 
approximation Q 

9 (pmo) Peoo Rm.s. X 2 

- 30° 26.92 31.47 320.42 6.0594 

30 -45° 14.81 17.86 476.45 5.0713 

45 - 60° 8.59 10.15 770.07 8.3343 

make juxtaposition with celestial coordinates as 
for functions ([[]) and (Q. 

III. RESULTS 

Revised parameters of individual lateral dis- 
tribution functions in Greisen's (([T]) and ([2])) ap- 
proximation showed no correlation neither with 
Galactic plane, nor with the region of UHECR 
region excess. From table [j] it is seen, that in- 



creased Ab in the region of interest (1.7h < a < 
3.7h) is not significant and the whole picture is 
spoiled by poor statistics. 
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Figure 4: Average scaling LDFs for three different 
zenith angle intervals. 

Difficulties in estimation of R m s did now 
allow us to use scaling approximation ([3]) in 
such analysis. In the work by MSU EAS 
group [5] authors have faced similar obstacles 
in R ms determination. It is worth mention- 
ing, that KASCADE-Grande group successfully 
used scaling formalism for estimation of muon 
density in air showers [6J. Besides, scintillation 
detectors used at the Yakutsk array may lead to 
sloping of charged particle distribution caused 
by registration of atmospheric muons and elec- 
trons from muon decay. If we consider this 
fact together with zenith- angular dependence of 
i? m . s . more closely, we can obtain more plausible 
estimation of this parameter. 
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